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Abstract 

We study the statistical properties of the scattering matrix associated with generic quantum graphs. The scat- 
tering matrix is the quantum analogue of the classical evolution operator on the graph. For the energy-averaged 
spectral form factor of the scattering matrix we have recently derived an exact combinatorial expression. It is 
based on a sum over families of periodic orbits which so far could only be performed in special graphs. Here we 
present a simple algorithm implementing this summation for any graph. Our results are in excellent agreement 
with direct numerical simulations for various graphs. Moreover we extend our previous notion of an ensemble 
of graphs by considering ensemble averages over random boundary conditions imposed at the vertices. We show 
numerically that the corresponding form factor follows the predictions of random-matrix theory when the num- 
ber of vertices is large — even when all bond lengths are degenerate. The corresponding combinatorial sum has 
a structure similar to the one obtained previously by performing an energy average under the assumption of 
incommensurate bond lengths. 
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1. Introduction 

Quantum graphs have recently attracted a lot 
of interest . A review paper containing a list 
of references to previous work can be found in 
|l) . The attention is due to the fact that quantum 
graphs can be viewed as typical and simple exam- 
ples for the large class of systems in which clas- 
sically chaotic dynamics implies universal spec- 
tral correlations in the semiclassical limit |p|,|l0|. 
Up to now we have only a very limited under- 
standing of the reasons for this universality. In a 
semiclassical approach to this problem the main 
stumbling block is the intricate interference be- 
tween the contributions of (exponentially many) 
periodic orbits [0,|2|. Using quantum graphs as 
model systems it is possible to pinpoint and iso- 
late this central problem. In graphs, an exact 
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trace formula exists which is based on the peri- 
odic orbits of a mixing classical dynamical sys- 
tem [0,|3|. Moreover the orbits can be specified 
by a finite symbolic code with Markovian gram- 
mar. Based on these simplifications it is possi- 
ble to rewrite the spectral form factor or other 
two-point correlation functions in terms of a com- 
binatorial problem This combinatorial prob- 
lem has been solved with promising results for se- 
lected graphs: It was shown that the form factor, 
ensemble averaged over graphs with a single non- 
trivial vertex and two attached bonds (2-Hydra) 
coincides exactly with the random-matrix result 
for the 2 X 2 CUE §. In §] the short-time ex- 
pansion of the form factor for A^-Hydra graphs 
(i. e. one central node with A^ bonds attached) 
was computed in the limit N —> oo, and in 0] 
a periodic-orbit sum was used to prove Anderson 
localization in an infinite chain graph with ran- 
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domized bond lengths. In g the form factor of bi- 
nary graphs was shown to approach the random- 
matrix prediction when the number of vertices 
increases. In the present paper we develop a gen- 
eral method to implement the combinatorial sum 
on a computer. Our results are always compared 
with direct numerical simulations. 

Most of the studies on quantum graphs pre- 
sented up to now were done with a fixed set of 
vertex boundary conditions at the vertices. In 
this case, one of the main premises to find spectra 
following the random-matrix predictions is that 
all bond lengths of the graph are rationally inde- 
pendent. In contrast, we will consider here also 
statistical properties of graph spectra averaged 
over the set of boundary conditions. We show 
that in this case the form factor is in agreement 
with the random-matrix theory (RMT) predic- 
tions even when all bond lengths are degenerate. 
Moreover, we show that one can express the form 
factor averaged over the vertex-scattering matri- 
ces as a combinatorial sum over families of orbits. 
This sum has the same structure as obtained in 
1^ by performing a spectral average. 

This paper is structured in the following way. 
In the following Section 2, the main definitions 
and properties of quantum graphs are given. We 
concentrate on the unitary bond-scattering ma- 
trix Sb which can be interpreted as a quantum 
evolution operator on the graph. Section 3 deals 
with the corresponding classical dynamical sys- 
tem. In Section 4, the statistical properties of 
the eigenphase spectrum of the bond-scattering 
matrix Sb are analyzed and related to the peri- 
odic orbits of the classical dynamics. Finally, our 
conclusions are summarized in Section 5. 

2. Quantum Graphs: Basic Facts 

We start with a presentation and discussion of 
the Schrodinger operator for graphs. Graphs con- 
sist of V vertices connected by B bonds. The va- 
lency Vi of a vertex i is the number of bonds meet- 
ing at that vertex. The graph is called v-regular 
if all the vertices have the same valency v. When 
the vertices i and j are connected, we denote the 
connecting bond by 6 = (*7 j)- The same bond can 
also be referred to as 5 = {Min{i, j), Max{i,i)) or 



b = {Max{i,j),Min{i,j)) whenever we need to 
assign a direction to the bond. A bond with coin- 
ciding endpoints is called a loop. Finally, a graph 
is called bipartite if the vertices can be divided 
into two disjoint groups such that any vertices 
belonging to the same group are not connected. 

Associated to every graph is its connectivity 
(adjacency) matrix Cij. It is a square matrix 
of size V whose matrix elements Cij are given in 
the following way 

„ _„ _ ( 1 if i,j are connected 1 
- - I otherwise J ' 

For graphs without loops the diagonal elements of 
C are zero. The connectivity matrix of connected 
graphs cannot be written as a block diagonal ma- 
trix. The valency of a vertex is given in terms 
of the connectivity matrix, by Vi = 
and the total number of undirected bonds is 

= 2 Sij = l 

For the quantum description we assign to each 
bond b — a coordinate Xi_j which indi- 

cates the position along the bond. Xij takes the 
value at the vertex i and the value Lij = Lj^i 
at the vertex j while Xj^i is zero at j and Lij 
at i. We have thus defined the length matrix 
Li^j with matrix elements different from zero, 
whenever Ctj ^ and L^j = Lj^i for b = 
1, B. The wave function ^> contains B compo- 
nents ^ffci (xfej, ^E* 62 (xfc^ ),..., ^'hg(xfcg) where the 
set {bi}f^i consists of B different undirected 
bonds. 

The Schrodinger operator (with h — 2m = 1) 
is defined on a graph in the following way: On 
each bond b, the component of the total wave 
function ^E" is a solution of the onc-dimensional 
equation 

("'^"^0 = fc2*6(a;) . (1) 

We included a "magnetic vector potential" Ab 
(with 5Re(Ab) ^ and = -^t) which breaks 
the time reversal symmetry. In most applications 
we shall assume that all the Af,'s are equal and 
the bond index will be dropped. On each of the 
bonds, the general solution of (|^) is a superposi- 
tion of two counter propagating waves 

= a,je'('=+-^''^)^''^ + aj- ,e'('=+^^-')^^-' (2) 
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The coefficients aij form a vector a = (oj^^ , 

• ■ •, ag^ , a'^_^ , • ■ ■ 1 '^7 complex numbers 
which uniquely determines an element in a 
2B— dimensional Hilbert space. This space corre- 
sponds to "free wave" solutions since wc did not 
yet impose any conditions which the solutions of 
(P have to satisfy at the vertices. 

The most general boundary conditions at the 
vertices are given in terms of unitary Vj x Vj 

vertex-scattering matrices cr;'""'2j(fc), where I and m 
go over all the vertices which are connected to j. 
At each vertex j, incoming and outgoing compo- 
nents of the wave function are related by 



^3.1 



kL-i 



(3) 



which implies current conservation. The particu- 
lar form 



2 

Vi 



(4) 



for the vertex-scattering matrices was shown in 
Q to be compatible with continuity of the wave 
function and current conservation at the vertices. 
(||) is referred to as Neumann boundary condi- 
tions. 

Stationary states of the graph satisfy (|^) at 
each vertex. These conditions can be combined 
into 



a = SB{k)a, 



(5) 



such that the secular equation determining the 
eigenenergies and the corresponding eigenfunc- 
tions of the graph is of the form |Q| 

det[/-S'i3(fc,A)] =0. (6) 

Here, the unitary bond-scattering matrix 

SB{k,A)^D{k-A)T (7) 

acting in the 2i?-dimensional space of directed 
bonds has been introduced. The matrices D and 
T are given by 



(8) 



T contains the topology of the graph and is equiv- 
alent to the complete set of vertex-scattering ma- 
trices, while D contains the metric information 
about the bonds. 

It is instructive to interpret the action of Sb on 
an arbitrary graph state \1/ as its time evolution 
over an interval corresponding to the mean bond 
length of the graph such that 



a(n) = S'ga(O), n = 0,1,2,, 



(9) 



Clearly the solutions of (H) are stationary with re- 
spect to this time evolution, n in represents a 
discrete (topological) time counting the collisions 
of the particle with vertices of the graph. 

3. Periodic orbits and classical dynamics 
on graphs 

In this section we discuss the classical dynam- 
ics corresponding to the quantum evolution (||) 
implied hy Sb- To introduce this dynamics we 
employ a Liouvillian approach, where a classical 
evolution operator assigns transition probabilities 
in a phase space of 2B directed bonds 0]. If pb{t) 
denotes the probability to occupy the (directed) 
bond b at the (discrete) topological time t, we can 
write down a Markovian Master equation of the 
form 



Pb{t + l) = ^UbM'pb'{t). 



(10) 



The classical (Frobenius-Perron) evolution oper- 
ator U has matrix elements 



tj.nm 



- s ■ p'-^^ 



(11) 



with Pljl^_^^j, denoting the transition probability 
between the directed bonds b = {j, i) and b' — 
{i,j')- To make the connection with the quantum 
description, we adopt the quantum transition 
probabilities, expressed as the absolute squares 
of matrix elements of Sb 



P, 



Note that PjSj' ^^'^ ^ ^^t involve any metric 
information on the graph. In general, they may 
depend on the wave number k. 



(12) 
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The unitarity of the bond-scattering matrix Sb 
guarantees Ub,b' = 1 and < UbM' < 1, 

so that the total probabihty that the particle is 
on any bond remains conserved during the evo- 
lution. The spectrum of U is restricted to the 
interior of the unit circle and i/i = 1 is always 
an eigenvalue with the corresponding eigenvector 
1) = {1,1, . In most cases, the eigen- 

value 1 is the only eigenvalue on the unit circle. 
Then, the evolution is ergodic since any initial 
density will evolve to the eigenvector |1) which 
corresponds to a uniform distribution (equilib- 
rium). The rate at which equilibrium is ap- 
proached is determined by the gap to the next 
largest eigenvalue. If this gap exists, the dynam- 
ics is also mixing. Graphs are one dimensional 
and the motion on the bonds is simple and stable. 
Ergodic (mixing) dynamics is generated because 
at each vertex a (Markovian) choice of one out 
of V directions is made. Thus, chaos on graphs 
originates from the multiple connectivity of the 
(otherwise linear) system 

Despite the probabilistic nature of the classi- 
cal dynamics, the concept of a classical orbit can 
be introduced. A classical orbit on a graph is an 
itinerary of successively connected directed bonds 
(*i)*2), (»2,*3), • ■ •■ An orbit is periodic with pe- 
riod n if for all k, {in+k, Wfc+i) = {ik,ik+i)- For 
graphs without loops or multiple bonds, the se- 
quence of vertices ii, Z2, • • • with i„i £ [1, V] and 
Ci^^i^^j = 1 for all m represents a unique code 
for the orbit. This is a finite coding which is gov- 
erned by a Markovian grammar provided by the 
connectivity matrix. In this sense, the symbolic 
dynamics on the graph is Bernoulli. This analogy 
is strengthened by further evidence: The number 
of n— PO's on the graph is ^trC", where C is the 
connectivity matrix. Since its largest eigenvalue 
Tc is bounded between the minimum and the 
maximum valency i.e. mhiVi < Fc < m.ax.Vi, pe- 
riodic orbits proliferate exponentially with topo- 
logical entropy « logFc. 

4. The spectral statistics of Sb 

We consider the matrix SB{k, A) defined in (M), 
(^). The spectrum consist of 2B points e'^'*"^^ 
confined to the unit circle (eigenphases). Uni- 



tary matrices of this type are frequently studied 
since they are the quantum analogues of classi- 
cal, area preserving maps. Their spectral fluc- 
tuations depend on the nature of the underlying 
classical dynamics The quantum analogues 
of classically integrable maps display Poissonian 
statistics while in the opposite case of classically 
chaotic maps, the eigenphase statistics conform 
with the results of RMT for Dyson's circular en- 
sembles. To describe the spectral fluctuations of 
Sb we consider the form factor 

K{n/2B)^^{\i,Sl\^) (n>0). (13) 

The average (. . .} will be specified below. RMT 
predicts that K{n, 2B) depends on the scaled 
time T = ^ only jl^, and explicit expressions 
for the orthogonal and the unitary circular en- 
sembles are known [ p^ . 

Using (j?]), (||) we expand the matrix products 
in tvSg and obtain a sum of the form 

trS'S(fc) = ^pe'^'^-^p+^'p) . (14) 

In this sum p runs over all closed trajectories on 
the graph which are compatible with the con- 
nectivity matrix and which have the topological 
length n, i. e. they visit exactly n vertices. For 
graphs, the concepts of closed trajectories and pe- 
riodic orbits coincide, hence ( p^ can also be in- 
terpreted as a periodic-orbit sum. From ( p^ it is 
clear that K{n/2B) = as long as n is smaller 
than the period of the shortest periodic orbit. 
The phase associated with an orbit is determined 
by its total (metric) length Lp = X]&ep and by 
the "magnetic flux" through the orbit. The latter 
is given in terms of its total directed length lp if we 
assume for simplicity that the magnitude of the 
magnetic vector potential is constant \Ab\ = A. 
The complex amplitude of the contribution from 
a periodic orbit by the product of all the elements 
of vertex-scattering matrices encountered 

■^-.-n<'Si.,„-n(<'i?)"''""'. (15) 

i. c. for fixed boundary conditions at the ver- 
tices it is completely specified by the frequencies 
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isometric orbits which can be used to write the 
result of fFel) in the form i,|6|j7l 



Fig. 1: Form factor of Sb for a double-loop graph (or 
equivalently Sg for a 4- Hydra). The full line shows ( |l8[ ) 
corresponding to zero magnetic field in the loops or two 
pairs of degenerate bond lengths in the Hyd ra, as shown in 
the inset. Dashed line (obtained from (bol)): A magnetic 
flux threads the loops and is averaged over. This corre- 
sponds to 4 incommensurate bond lengths of the Hydra 
graph. 



ripir, s, t) of all transitions (r, s) — * (s, t) . Insert- 
ing (|l^) into the definition of the form factor we 
obtain a double sum over periodic orbits 

Kin/2B) = ^ ( E -^P^l' ^ 

exp {ifc(Lp -Lp,) + iA{lp - Ip,)} ^ . (16) 

This will be our starting point for the combina- 
torial approach presented in the following subsec- 
tion. 

4.1. Energy Average 

In previous work the average in ( p^ was 

always taken with respect to the wave number k 
and, if present, also over the magnetic vector po- 
tential A. Within the present subsection we will 
stick to this procedure. Provided that the bond 
lengths of the graph are rationally independent 
and that a sufficiently large interval is used for 
averaging, only terms with Lp — Lp> and Ip = Ipi 
survive. A completely equivalent result can be 
obtained by regarding the bond lengths Lb of the 
graph as random numbers and performing an en- 
semble average over them. 

Note that Lp = Lpi does not necessarily imply 
p — p' or that p,p' are related by some symme- 
try because there exist families C of distinct but 



K{n/2B) = 



(17) 



The outer sum is over the set Tn of families, while 
the inner one is a coherent sum over the orbits be- 
longing to a given family (= metric length), ( p^ ) 
is exact, and it represents a combinatorial prob- 
lem since it does not depend anymore on metric 
information about the graph (the bond lengths). 

Let us illustrate the application of (|l^ in 
the following simple system: We consider a Hy- 
dra graph with four arms and impose Neumann 
boundary conditions at the central node such 
that forward and backward scattering amplitudes 
are given according to (||) by cr/ = 1/2 and 
(T(, = —1/2. For the moment we require that there 
are two pairs of identical bond lengths La, Lb 
(see Fig. 1). Since Hydra is a bipartite graph, 
an orbit can return only after an even number of 
vertex-scattering events (note that the dead ends 
of the arms are considered as vertices). It is then 
convenient to replace the scattering matrix in ( [T^ ) 
by S^. As one can see from (|l^, that the form 
factor obtained in this way is equivalent to the 
one obtained from a graph which is formed by two 
loops of lengths 2La and 2L b connected by a sin- 
gle vertex (inset of Fig. 1). Consider now the peri- 
odic orbits contributing to K{n/2B). The length 
of such an orbit is given by L — ^iaLa + ubLb 
where ua, tlb denotes the number of traversals 
of A and B loop, respectively, and n = ua + nB- 
If La, Lb are rationally independent, a family of 
isometric orbits in (|l^ contains all orbits shar- 
ing UA, riB- These orbits differ, however, in their 
symbolic code by the order of A's and B's and by 
the orientation in which the loops are traversed. 
This is reflected in the amplitudes Ap of the or- 
bits: If two consecutive loops A A have different 
orientation, the amplitude Ap contains a factor 
(Tb, while the same orientation results in cr/. Dif- 
ferent loops following each other (AB or BA) al- 
ways result in a factor af. Suppose now there is 
an orbit which contains in its code a sequence of 
the form • • AB1B2B5 ■■ - A- ■■. Then there will be 
another orbit ■ • ABiB'2B's ■ ■ ■ A - ■ ■ which differs 
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from the first one by an inversion of the orienta- 
tion (denoted by the prime) of all but the first 
loop B. It is easy to see that the corresponding 
amplitudes cancel. This results in a tremendous 
simplification of the problem: It suffices in ( |l7| ) 
to keep track of all trajectories for which the can- 
cellation does not apply. These are of the form 
A", B", (AB)"/2 or (BA)"/^ the latter two ex- 
isting for even n only. Each orbit A" contributes 
the amplitude 1/2" irrespective of the orienta- 
tions of the n loops. On the other hand there 
are 2" possibilities to prescribe these orientations 
and hence the total contribution from each of the 
families A" and B" is 1. There are also 2" trajec- 
tories of the form (AB)"/^ and (BA)"/^, respec- 
tively, each with amplitude 1/2". However these 
orbits belong to the same family jia ~ ns = n/2, 
and their amplitudes must be added coherently 
in (17). This results in the total contribution 4 



from this family. After adding the contributions 
from the three different families and normalizing 
with 2B = 4 we are left with 



K{n/2B) = 



1/2 
3/2 



1,3, 
2,4, 



(18) 



which is shown in Fig. 1 with a solid line. Ne- 
glecting the point n = 0, the form factor sim- 
ply oscillates around the value 1 — even within the 
Heisenberg time r = n/2B = 1. This is reminis- 
cent of the behavior of a Poissonian spectrum. 
The reason for this lack of correlations is that 
there exist regular sequences of eigenstates which 
are completely confined to one of the two loops or 
one pair of degenerate bonds, respectively. These 
states disappear when all arms of the Hydra are 
incommensurate, and the corresponding form fac- 
tor (dashed line in Fig. 1) clearly reflects this 
change in the properties of the spectrum. 

The arguments which led to (|l^) are deceiv- 
ingly simple. In general, the combinatorial prob- 
lem (I?!) is very hard and cannot be solved in 
closed form. Even for the 4-Hydra with incom- 
mensurate bond lengths this seems to be the case. 
Nevertheless exact result for finite n as those 
shown in Fig. 1 with a dashed line can always 
be obtained from (^^ using a computer algebra 
system such as Maple |l^. This will be shown in 
the following. To prescribe topology and bound- 
ary conditions, one should provide the time evo- 



lution operator Sb according to (|7|) but with the 
phases e''^^' J left as unspecified variables. For 
example, the operator Sb for the graph discussed 
above has the form 



1/2 <t,B 
1/2 4>B 
1/2 4>B 
-1/2 4>B 



will be as many phases as there are incommen- 
surate lengths in the graph. The quantity trS'g 
can now be represented as a multivariate polyno- 
mial of degree n in the variables 0^, i. e. 





' -1/2 


1/2 4) A 


1/2 4>B 




1/2 


-1/2 4>A 


1/2 4>B 




1/2 (Pa 


1/2 4, A 


-1/2 ct>B 




1/2 4> A 


1/2 4, A 


1/2 4>B 


where (I)a/b — 




In the g( 



tr^S 



cv 



(19) 



where Vn runs over all partitions of n into non- 
negative integers n = pi + p2 + ■ ■ ■■ The form 
factor is then simply given as 



Kin/2B) - J2 



(20) 



The task of finding the coefRcients c-p can be ex- 
pressed in Maple with standard functions, such 
that after the initialization of Sb the following 
four simple lines represent a completely general 
algorithm for the exact computation of the form 
factor of an arbitrary graph: 

sn:=evalm(S~n) ; 

tn : =expand (trace ( sn) ) : 

c: = [coeff s(tn)] ; 

K : =suin(c [p] "2 ,p=l . .nops(c) ) ; 

In practice, however, the computation is re- 
stricted to the first few n since the numerical ef- 
fort grows exponentially fast. In Fig. 2 we com- 
pare the results of ( |20| ) with direct numerical av- 
erages for regular (fully connected) graphs with 
V = A and V — h vertices with and without mag- 
netic field breaking the time-reversal symmetry. 
The results agree indeed to a high precision. Al- 
though this could be regarded merely as an ad- 
ditional confirmation of the numerical procedures 
used in |l| , we sec the main merit of (^0|) in being 
a very useful tool for trying to find the solution 
of (|l7f) in closed form. 
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FIG. 2. Form factor of Sb for regular graphs with 
V = A vertices (top) and V = 5 vertices (bottom). In the 
right panels an additional magnetic field destroyed time- 
reversal symmetry. Circles: exact results obtained from 
(hj). Stars: numerical average over 2,000 values of fc. The 
solid line is to guide the eye. The prediction of the ap- 
propriate random matrix ensemble are shown with dashed 
lines. 



4.2. Ensemble Average over the Boundary 
Conditions at Vertices 

In the present section we generate an ensemble 
of graphs by randomizing the vertex-scattering 
matrices a^^\ The length matrix L which con- 
tains the lengths of the bonds and the connec- 
tivity matrix (topology of the graph) are kept 
constant. Moreover all the bond lengths are 
equal Li, = 1. We computed the form factor 
for fully connected graphs. Our results for val- 
ues of y = 10, 15, 20 are presented in Fig. 3. The 
RMT two-point form factor ||T^ is also displayed 
in Fig. 3 for comparison. The results show quite 
a good agreement with the predictions of RMT 
for the circular ensembles. Moreover, one can see 
that as the number of vertices V increases, the 
agreement with RMT becomes progressively bet- 
ter. The deviations from the smooth curves are 
not statistical, and cannot be ironed out by fur- 
ther averaging. Rather, they are due to the fact 
that the graph is a dynamical system which can- 
not be described by RMT in all detail. 

In the remainder of this subsection, we would 
like to point out the possibility for turning (16) 
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Fig. 3: The form factor of the eigenphase spectrum of 
Sb for fully connected graphs with V = 10, 15, 20. Bold 
dashed lines are the expectations for the COE and CUE 
expressions. The data are averaged over realizations of 
the (t's, as explained in the text. All the bond lengths are 
kept constant and equal to j = 1. 



into a combinatorial problem using the ensemble- 
average over tr'^*^ . We assume that cr'-'-' i = \, . . .V 
are independent random matrices chosen from the 
CUE. In this case, the averaging in ( jl6| ) factorizes 
with respect to the vertices and we find 



s [r,t] 



(21) 



Averages over the CUE of this type have been 
computed in the literature, see e. g. iQ. The 
result is of the form 



n 



p,p' 



(22) 



where P, P' run over the permutations of 1, . . . , n 
and Cp-ipi are some constants which can be 
obtained from a recursion relation. Applied to 
( |2l] ) this implies, that the contribution from a 
pair of orbits p,p' vanishes unless for all s,t 

rip/ (r, s, t), i. e. unless the frequencies of tra- 
versal of any directed bond s — > t coincide for p 
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and p' . For the case of a graph with time-reversal 
symmetry, i. e. when the vertex-scattering matri- 
ces are taken at random from the COE, a sim- 
ilar result ensures that only those pairs of or- 
bits survive the averaging in (|l]) which agree 
in the traversals of all undirected bonds. This 
means that the family structure underlying the 
set of periodic orbits is exactly the same as in 
the case of a graph with incommensurate bond 
lengths but fixed boundary conditions: The fam- 
ilies are formed by those orbits which differ only 
in the time order of the traversed bonds. Sub- 
stituting this result into (|l^) we see that all os- 
cillating phases containing metric information on 
the graph drop independent of the precise values 
of the bond lengths and we are indeed left with a 
combinatorial problem. 



5. Conclusions 

In this paper we have tried to contribute to 
the understanding of the statistical properties of 
the unitary quantum time evolution operator de- 
rived from quantum graphs. This problem is rel- 
evant since it is a paradigm for the as yet unan- 
swered question precisely under which conditions 
the quantum analogues of classically chaotic sys- 
tems are universal and follow the random-matrix 
predictions. Fully connected quantum graphs 
show this universality when the number of bonds 
becomes large. In extension to previous work we 
have demonstrated this result in the case where 
an ensemble of graphs is introduced by randomiz- 
ing the boundary conditions at the vertices. The 
corresponding ensemble average can replace the 
previously considered spectral average, and the 
RMT results are in this case approached even 
when all the bond lengths are incommensurate. 
We have the hope that this kind of ensemble av- 
erage might provide an easier access to an analyti- 
cal treatment of the spectral properties of graphs. 
One possible approach of this goal is the use of 
combinatorial methods to perform the periodic- 
orbit sums related to spectral two-point correla- 
tions. 
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